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1. Introduction 

The infrared behavior of propagators has received a lot of interest in recent years. Many of the 
discussions have evolved around the zero momentum value of the gluon propagator and the en- 
hancement of the ghost. Until recently, the lattice calculations and analytical approaches agreed 
on a vanishing gluon propagator at zero momentum and a ghost propagator which is infrared en- 
hanced. One particular model, the Gribov-Zwanziger framework, which takes into account Gribov 
copies, also agreed on this behavior. However, in the last couple of years, the understanding of 
the behavior of the propagators has changed dramatically. The common believe is now that the 
ghost propagator is not enhanced and the gluon propagator attains a finite nonvanishing value at 
zero momentum. Such propagators were found by lattice as well as by many analytical approaches 
[jl|, Q |5[ Qj. In particular, in the Gribov-Zwanziger framework, the new behavior of the ghost 
and gluon propagator was explained by taking into account the existence of a certain dimension 2 
condensate. This framework was called the refined Gribov-Zwanziger framework [Q, 

The condensate of interest in the refined GZ action stems from a BRST invariant operator of di- 
mension d = 2. However, as is well known, the GZ action breaks BRST symmetry, so there is no 
reason to a priori investigate only BRST invariant operators. In fact, more d = 2 condensates can 
be investigated without spoiling the renormalizability. These condensates are dynamical objects, 
giving rise to dynamical mass scales. To get an estimate for the latter, we have to compose the 
effective potential and look for the minimum of this potential. We shall do this for the first time, 
our aim being that of providing evidence that these condensates are in fact present. 

2. The Gribov-Zwanziger action 

Let us start by reviewing the Gribov-Zwanziger action. The Yang-Mills action is given by 

Sym = ^Jd 4 xF£ v F« v . (2.1) 

However, the Yang-Mills action needs to be gauge fixed in order to define the path integral, as 
else, it is ill-defined. This can be done with the Faddeev-Popov action which, in the Landau gauge, 
d^A^ = 0, reads 

Sfp = Sym + , (2.2) 

with 

% = J d 4 x (><V£ +rd^Dfc h ) , (2.3) 
being the gauge fixing part. This action enjoys BRST symmetry, given by 

sA a pi = -(D li c) a , sc a = X -gf abc c b c c , sc" = b a , sb a = 0, s 2 = 0, (2.4) 

which is important to prove perturbative unitarity and renormalizability. However, as shown by 
Gribov, this gauge fixing suffers from Gribov copies [§]. If we take e.g. the Landau gauge, d^A^ = 
0, it is very easy to show that there exist gauge equivalent fields A'^ also fulfilling the Landau gauge 
d^A' = when having zero modes of the Faddeev-Popov operator 

J?*> = -d^Df = -d^d^-gf^) . (2.5) 
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Therefore, Gribov proposed to restrict the region of integration in the path integral to Q, 

Q={A jU |a M A M = 0,^>0)} , (2.6) 

the region which is free of zero modes of the Faddeev-Popov operator ^ ah . Unfortunately, it is 
not so difficult to show that there are still Gribov copies inside the Gribov region, as first proven 
in [10]. In 1989, Zwanziger implemented this restriction to the Gribov region to all orders: the 
Gribov-Zwanziger action was born [11]. We present immediately the localized version of this 
action, namely 

Sqz = Spp + So + S r , (2-7) 
with Spp the usual Faddeev-Popov action, 

S FP = I J d 4 xF« v F; v + J d 4 x (V<V£ +rd^Dfc b ) , (2.8) 

and So and S r given by 

So = J d A x (^fJt"* 3 ^ - o^Jt ab (i>f) , 

Sy = -fg J d 4 x (f^itf + qftW + ~ g {N 2 - 1) T 2 ) • 

The fields (<p^ c , <p^ c ) are a pair of complex conjugate bosonic fields, while (f0^ c , G)™) are anticom- 
muting fields. 7 is not free, but fixed by the horizon condition, 

dr 

whereby T is the quantum action defined as, 

e- r = j[d<l>]e- s , (2.10) 

and [dtp] represents the integration over all the fields. In the case of the local action, this condition 
is equivalent to, 

{gfckAi « = 2y 2 4(n 2 - 1) . (2.1 1) 

Looking at the BRST transformation of the new fields 1 , 

5<pf = 0?, sft)f = 0, scSf =7p° , syf = 0, (2.12) 
one can check very easily that Sqz is no longer invariant under the BRST symmetry s, 



sS GZ = s (S YM + S gf + So + S 7 )=s (S r ) =gfj d d xf bc (a£o)J c - (D™c m ) (V* c + % 



(2.13) 



Therefore, this action breaks the famous BRST symmetry which, among other things, leaves open 
the question of unitarity. For more on this breaking, we refer to [0]. Let us mention that when re- 
moving the Gribov horizon, i.e. setting 7 = 0, there is no breaking and the extra fields (<p,lp, ft), ft)) 



1 We shall use a shorthand notation for the new fields: e.g. (p^f = (p" 
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give rise to a BRST quartet and decouple from the theory. It is remarkable that the breaking of the 
BRST does not spoil the renormalizability of the GZ action, due to the soft nature of the breaking. 

Let us have a look at the ghost and gluon propagators predicted by the Gribov-Zwanziger model 
[§, |ll|, 12, 13p. Firstly, the gluon propagator is positivity violating, infrared suppressed and vanish- 



ing at zero momentum. E.g. the tree level gluon propagator is given by, 



(Al(p)K(-p)) = S* (V - P -f-) p4 + P " g2Nf , (2.14) 



whereby also at higher loops, the gluon propagator becomes zero at zero momentum [13]. Sec- 
ondly, the ghost propagator is enhanced at zero momentum, 

limG(£ 2 )~A, (2-15) 

which has also been explicitly checked up to two loops, [12, 13]. However, both results do not 
seem to agree with the lattice results, see [[l], 

3. Refining the Gribov-Zwanziger action 



As the horizon condition ( |2- 1 1[ ) is equivalent with a dimension two condensate, it would appear 
naturally that other such condensates could appear in the Gribov-Zwanziger approach. In [Q], we 
have investigated the following two condensates, (A 2 ^ and (q>"(p" — Wfco"), which gave rise to the 
refined Gribov-Zwanziger action, 

Srgz = S cz + J d 4 * Qm^JAJ - M 2 (<p?<p? - . (3.1) 

The masses m 2 and M 2 associated with the condensates have a dynamical origin. 

The propagators of the refined action display a different behavior now. Firstly, the gluon prop- 
agator remains infrared suppressed and positivity violating, but does not vanish anymore at zero 
momentum. At tree level it is given by 

D( 2) = p 2 +M 2 

{P ' p* + (M 2 + m 2 )p 2 + 2g 2 Nf+M 2 m 2 ' K ' 

and it remains non zero at higher loop order Secondly, the ghost propagator behaves oc \/k 
and is thus not enhanced anymore. The behavior of the propagators is in qualitative agreement with 
the lattice data of [[I], |14| ]. For some quantitative numbers, we refer to [^], whereby the tree level 
gluon propagator Q3.2[) was fitted to the lattice data. The following estimates were found, 



,2x2' 



m 2 p 


a -1.78 GeV 2 


M 2 p 


a 2.14 GeV 2 , 


A 4 p 


a 0.26 GeV 4 


D(0) p 


a 8.2 GeV 2 , 


Jl=10 GeV P 


a 3 GeV 2 , 



^) u =ioGev«3GeV^, (3.3) 
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in association with a very good quality of the fit, and with the value of (g 2 A 2 ) in the ballpark of 
other estimates from different approaches [ |l6[ f% . 

4. A further refining of the GZ action 



So far, the main ingredient in the refining of the GZ action was the inclusion of a dimension 
2 condensate for which the corresponding operator is BRST invariant, namely <p?<pf — co"co" = 
s(cof <pf). However, as we have explained, the GZ action breaks the BRST symmetry, so there is 
no reason to only investigate this d = 2 BRST invariant operator. In fact, not only are there more 
possible condensates which can be present, we shall also split the operator <pf<pf — ft)" ft)" into two 
separate pieces. All possible renormalizable d = 2 operators G{ in the GZ action, which have ghost 
number zero, are given by 

0i = {A^^^m^Wi^t} ■ (4.i) 

We choose to only investigate operators for which both color and Lorentz indices are fully con- 



tracted, e.g. like (p-'lp" = (p^Tp^ . However, as one can find in [|15[], there are other possibilities to 
combine the color indices. If one wants to be absolutely complete, one would have to take into 
account all possible color contractions. Unfortunately, doing this would be hopelessly complicated 
and we hope that we have captured the mayor physics by taking only one color combination. How- 
ever, in principle, different color combinations are possible. 



We have added all the d = 2 operators 0{ to the GZ action, and we have proven that this can 



be done in a renormalizable fashion [|18|]. As there is no lowest order coupling of ft) and ft) to 
the gluon sector, for the calculation of the effective action, we shall not take this condensate into 
account. The enlarged action is then given by 



-CGZ 



S G z + 



d 4 x 



2^ 



aQQ-xQz 



d 4 x 



^Gyflpt + ^Gtftf+pGG 



(4.2) 



whereby Q, X, G and G are the new sources coupled to the relevant operators, while £, a, % an d 
p are the dimensionless LCO parameters of the quadratic terms in the sources which are needed to 
account for the divergences present in the correlation functions like (0j(k)0j(—k)), with 0i one of 



the operators given in expression (6.1 ), see [19]. This action is multiplicative renormalizable, how- 
ever, the renormalization is non-trivial. Due to mixing of the sources Q and X, the renormalization 



requires some effort J18Q. 



Let us again discuss the propagators. The gluon propagator becomes slightly more complicated, 
although remains infrared suppressed and non-zero at zero momentum. The tree level expression 
is given by 



D(P 2 ) 



(M 2 + p 2 ) 2 -pp 



t 



MV + p 6 +M 2 (2p 4 + A 4 ) - A 4 ^^ +m 2 ( (M 2 + p 2 ) 2 - ppt ) + p 2 (X 4 - pp 



(4.3) 
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Notice that this propagator reduces to the form of the RGZ propagator ( fO| ) in the case that p = p 

M 2 + p 2 + p 



D(P Z ) 



p 4 +M 2 p 2 + p 2 {p+m 2 )+m 2 (M 2 +p) + X 4 ' 



(4.4) 



Comparison with existing lattice data shall determine whether this is the case or not. Secondly, the 
ghost propagator remains not enhanced, similar to the previous section. 

5. The effective action 

Fortunately, we shall be able to calculate the effective action at lowest order. For this, we shall 
follow the local composite operator (LCO) formalism developed in [|19|, ^]. In fact, it is interesting 
to notice that before is was not possible to apply the LCO formalism due to the specific form of the 
operator lp"(p" — ft)" ft)", as we were not able to perform a Hubbard Stratonovich transformation to 
find the effective action 2 . However, by splitting this operator as we did, it is possible to apply the 
LCO formalism. After a long analysis, we were able to calculate the effective action|] 



rd) = (^11 



(M 2 



16k 2 
-2{M 2 + pp t 

1 48 {N 2 - 
+ 2 53N 



PP 



t) 2 ln 



m 2 -7p7i 



w 



+ (M 2 + Jppt) 2 ln 



M 2 + ^/ppl 



w 



+ 



3(N 2 



64 n 2 

p 



i)r 5 ( 4 

- - (m 



2A 4 )+ 2 ln (zZil 2 ln (zZ?J 
P P 



P- 



l) 2 



2{N 2 
53 pi 



1) 



A 4 3 X A 



+ - 



Ng 2 2 32n 2 



(N 2 



V, 



+ 



9 N 2 - 1 m 4 
13 iV V 



I -Nit 

M 24 (N 2 - l) 2 J g 2 

24(N 2 -l) 2 M 4 161 N 2 



PP' 



Ira 4 



35 



N 



52 167T 2 2 



-M 4 0Ci +M 2 m 2 %i 



(5.1) 



whereby y\, yi and are the solutions of the equation y 3 + (m 2 + 2M 2 )y 2 + (A 4 +M 4 — pp f + 
2M 2 m 2 ).y+M 2 A 4 + l/2(p+pt)A 4 +M 4 »i 2 



m 2 pp f 



: and j4 and y$ of the equation y + 2M y + 
M 4 — pp T = 0. It is not only necessary to calculate all the one loop diagrams of the theory, also 
some two loop information is needed in order to determine the one loop effective action fll9| ^CQ . 
However, this remains for future work, so we have parameterized the unknown information into 
the three constants p\ , a\ and Xi ■ 



In order for the condensates to be indeed present, we have to show that the minimum of the ef- 
fective action requires a non zero value at least for one of the condensates (<p<p), ((pep) or ((pip). In 
order to show this, let us simplify the effective action by setting p = p^ = 0. This is equivalent with 
not considering the condensates {(pep) and {(pep) and only taking {(pip) into account. We remind 
the reader that already a non zero value of the condensate {(pip) is already sufficient for D(0) ^ 0. 
In this case, the effective action becomes: 



r (i) 



(W 2 



167T 2 



i) 2 r a m 2 

' 2M 4 ln =T 



-2M 1 



+ 



3(;v 2 -i; 

647T 2 



i( m 4_2A 4 )+M 4 ln^ 
6 p 



2 More details on this will be provided in JT^]. 
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+ylln ££> +y * ln t£l -artn^l - _,)*; + 



M M M 



Ng 2 2 327T 2 



9iV 2 -lm 4 24(/V 2 -l) 2 M 4 161^V 2 -lm 4 4 , , 

whereby y 2 and y 3 are given by ± ^-m 2 -M 2 ± Vm 4 - 2M 2 m 2 +M 4 - 4A 4 ' 
In order to find the minimum, we should derive this action w.r.t. m 2 and M z and put the equa- 



tions equal to zero. In addition, we should also impose the horizon condition (|2.9|). Therefore, we 
have the following three conditions, 

dr „ dr n dr 

0, ^=0, ^rr = 0, (5.3) 



dM 2 ' dm 2 ' dX 

which have to be solved forM 2 , m 2 and A 4 . Unfortunately, it is impossible to solve these equations 
exactly due to the two unknown parameters a\ and Xi ■ However, we can uncover if M 2 = can be 
a solution of the above expression. Looking at the three conditions, it is clear that we will obtain 
three equations with only two unknown parameters (m 2 , A 2 ). Although one of these equations still 
contains the unknown constant parameter %\, it is highly unlikely that this system of equations 
will have a solution for (m 2 , A 2 ). We take this as an indication that the condensate (<p<p) is indeed 
present, thereby suggesting the dynamical transformation of GZ into a refined GZ. 



6. Conclusion 



In this paper, we have made the refining of the GZ action more complete. We have shown that more 
d = 2 condensates are possible than the ones so far considered, i.e. (A^A^) and (<p"<pf — fflftof ). 
Moreover, as the GZ action breaks the BRST symmetry anyway, there is no reason to consider the 
operator Tp" q>" -co" cof as a whole. In summary, we have taken the following operators into account 
to further refine the GZ action: 



Gi = {A^,<pf<pf !<P f^,9f9f} , (6.1) 
while not spoiling the renormalizability. 



A second aim was to calculate the one loop effective action including all the condensates. We 
have succeeded in finding this effective action, see expression (5.1), although there are still some 
unknown constants stemming from two loop diagrams. Though, we are able to show that it is very 
likely that the minimum of the effective potential requires a non zero value of the condensate <pf<pf . 
This would prove the existence of a refined GZ action from a dynamical point of view. 
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